Kicked rotor quantum resonances in position space 
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We present an approach of the kicked rotor quantum resonances in position-space, based on its 
analogy with the optical Talbot effect. This approach leads to a very simple picture of the physical 
mechanism underlying the dynamics and to analytical expressions for relevant physical quantities, 
such as mean momentum or kinetic energy. The ballistic behavior, which is closely associated to 
quantum resonances, is analyzed and shown to emerge from a coherent adding of successive kicks 
applied to the rotor thanks to a periodic reconstruction of the spatial wavepacket. 



I. INTRODUCTION 

The kicked rotor is a simple system that plays a central 
role in studies of classical and "quantum chaos". The lat- 
ter is defined as the quantum behavior of a system whose 
classical counterpart is chaotic. In its simpler form, a 
kicked rotor (KR) is formed by a particle orbiting a fixed 
circular orbit to which an instantaneous force (a kick) 
is applied periodically. The corresponding classical dy- 
namics is found to be regular (periodic) for weak kick 
intensities, for intermediate kick intensities chaotic re- 
gions develop in limited zones of the phase space, and for 
strong enough forcing an ergodic diffusion appears [l| . In 
1995, Raizen and co-workers established this system 
as a privileged ground for studies of quantum chaos by 
realizing experimentally a kicked rotor in the quantum 
regime with laser-cooled atoms placed in a kicked (i.e. 
rapidly turned on and off) laser standing wave. 

Despite its apparent simplicity, the quantum kicked 
rotor (QKR) has very remarkable dynamical properties. 
One of the most studied is the so-called "dynamical lo- 
calization": in contrast to the classical case, the ergodic 
diffusion does not last forever in the quantum case. Af- 
ter a characteristic "localization time", the diffusion is 
stopped by destructive quantum interferences 0, [1, 
Another feature of the QKR that has been the object of 
a recent burst of theoretical and experimental activity is 
the existence of quantum resonances (QRs), whose most 
dramatic manifestation is the appearance, for specific val- 
ues of the parameters, of a ballistic motion, instead of a 
diffusive or a localized behavior. Quantum resonances 
are the main subject of the present paper. 

The experimental realization of the QKR has triggered 
in the last decade an impressive number of studies involv- 
ing dynamical localization 0, 0, 0, S, S lll|. q uantum 
transport [H, [H, [ll; [M . ratchets [lllllliil, chaos- 
20|, classical and quantum reso- 
Quantum resonances have been 
used in studies of fundamental aspects of quantum chaos 
such as quantum stabilization [25!] or measurements of 
the gravitation [2^. "High-order" quantum resonances 
were also observed recently [13, HE HI] both with a Bose- 
Einstein condensate and laser-cooled atoms. 

In the present work, we show that the kicked rotor 
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quantum resonances have a very simple and intuitive in- 
terpretation in position space, as opposed to the more 
common momentum space representation (see [sol ] and 
references therein). QRs in position space have been pre- 
viously considered by Izrailcv and Shepelyansky [J, Isij . 
Here, we extend their approach and build a simple phys- 
ical picture that enlightens the underlying physics and 
allows the calculation of quantities of experimental rele- 
vance, such as average momentum or the kinetic energy. 



II. PHYSICAL ORIGIN OF QUANTUM 
RESONANCES 



In this section we describe the fundamental properties 
of the kicked rotor dynamics and discuss the physical ori- 
gin of quantum resonances. We aim at a simple presen- 
tation that puts into evidence the physical mechanisms 
even for the reader unfamiliar with quantum resonances. 
Some aspects of the following discussion may be consid- 
ered as trivial, but the discussion is nevertheless neces- 
sary to introduce ideas and notation non-ambiguously. 

The atom-optics realization of a (quantum) kicked ro- 
tor consists in placing laser-cooled atoms in a far-detuned 
laser standing wave. In such conditions, the atoms feel 
the light intensity mechanical potential affecting 

their ccntcr-of mass degree of freedom [s^ [H, [13] . The 
standing wave is pulsed periodically, being on for a time 
interval so short that the motion of the atoms can be ne- 
glected; in such conditions, the pulses can be considered 
as delta functions (kicks). The standing wave sinusoidal 
modulation of intensity generates a spatial potential in 
sin(2fcix), where fei, = 27r/AL is the wavenumber of the 
radiation forming the standing wave. As the de Broglie 
wavelength of laser-cooled atoms is about A^/S, it is com- 
parable to the periodicity of the potential, Xl/2, and the 
system is in the quantum regime, provided that decoher- 
ence is negligible during the experiment. In practice, this 
implies using a far-detuned radiation to reduce sponta- 
neous emission to acceptable levels. 

We shall use a normalized spatial coordinate X = 
2kLX = 2x/Xl which plays in fact the role of a cyclic 
variable: the spatial periodicity of the potential implies 
that the physics is the same if translations by a multiple 
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of Ai,/2 are performed. We can thus use either a "hnear" 
("unfolded") representation of the KR, using the Hnear 
variable X or a "cyclic" or "folded" representation using 
the angular variable 6 = X[inod 27r] (a: [mod Ai/2] in 
usual units) fssj . 

Choosing, moreover, units such that the mass of the 
particle is 1 and the time-period of the forcing is T = 1, 
the Hamiltonian of the KR reads 



JV-l 



p2 

F = — + KcosX ^ 5{t - n) 

n=0 



(1) 



where P is the momentum (scaled by a factor M/2ki,T ) 
X the position in the periodic potential, K (usually called 
"stochasticity parameter") the intensity of the kicks, and 
n a discrete time corresponding to the rt*'* kick {n is an 
integer in normalized units). 

Labeling Xt,Pt the position and momentum imme- 
diately after the t*'' kick, and integrating the classical 
Hamilton equations of motion corresponding to Ec^. P]) 
produces the so-called Chirikov's "standard map" 

Xt+i = Xt + Pt 

Pt+i = Pt+KsinXt+i. (2) 
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The classical dynamics is found to be periodic below a 
critical value of the kick intensity Kc ~ 0.9716. Chaotic 
regions appear in phase-space for K > Kc, and progres- 
sively grow as K increases. For A' > 5 the islands of 
stable dynamics arc barely visible, and the classical dy- 
namics becomes an ergodic diffusion in phase space. The 
average kinetic energy then increases linearly with time 
(or kick number t): (P^) /2 = Dct, where « 
is a diffusion coefficient that can be explicitly calculated 

The standard map presents well-known classical reso- 
nances, also called accelerator modes. For example, set 
K = 27r, and consider a particle with initial conditions 
-''^0 = 7''/2 and Pq = 0. Iteration of Eqs. ^ shows that 
Pt = 2TTt and Xt = 2Ti{t^ - t) + n/2. The momentum 
increases linearly with time and the kinetic energy P^/2 
increases quadratically with time, which is a signature of 
a ballistic dynamics, in contrast with the linear increase 
in the diffusive case. The origin of the ballisticity is easily 
seen: for this particular choice of K and of the initial con- 
ditions, the particle is always kicked at the same position 
(modulus 27r), that is sinX(= sin(27r(t-^ — t) +tt/2)= 1, 
and thus receives the same amount of linear momentum 
per kick. For arbitrary values of K or of the initial con- 
ditions, the particle is kicked in different positions, and 
the effect of some kicks compensate the effect of other 
kicks, leading to a slower increase of the energy. One 
can easily convince oneself that, in general, there is a 
classical resonance for K = 2'Kp with integer p and ad- 
equate initial conditions. The classical resonance was 
experimentally observed in the atom-optics realization of 
the KR [2]| . There are also classical antiresonances: e.g. 
for the above initial condition and K = 3tt/2, successive 



kicks have opposite directions, and the momentum jumps 
endlessly between two values, Pt = and Pt+i = 37r/2. 

In order to study the quantum dynamics of the KR, one 
considers the one-period evolution operator, also called 
Floquet operator: 
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where k = Ahk\T/M is the normalized Planck's constant 
resulting from the definition of normalized variables sat- 
isfying the commutation relation [X, P] = ik; it thus 
describes the "quanticity" of the system {k^ is the 
classical limit). The above operator relates the quantum 
state after the t*'* kick to the quantum state after the 
{t— 1)*'' kick [33|. It is a remarkable fact that this opera- 
tor factorizes into the product of two exponentials, a free 
evolution followed by the kick effect, despite the fact that 
X and P do not commute. This is a consequence of the 
(5-function time dependence: during the instantaneous 
kick, the evolution related to the kinetic energy term is 
negligible. If one starts e.g. in the P-representation, 
the free evolution (from t+ to {t + corresponds to 
simply adding a phase. One can then convert to the X- 
representation where the kick operator is diagonal and 
also simply adds a phase; one then goes back to the P- 
representation. The calculation of the quantum evolution 
over a period thus "costs" only two Fourier transforms 
and two multiplications. It is this formal and numerical 
simplicity that made the QKR so popular. Despite this 
simplicity, the cross-action of these two operators makes 
the QKR dynamics very rich: the free evolution operator 
mixes space components and the kick operator adds new 
momentum components, generating a complex quantum- 
interference pattern. 

An important property of the kick operator is that it 
is periodic in space. This means that its cigenstates have 
a Bloch-wave (BW) structure. This is formally seen by 
developing it in terms of Bessel functions (noting that 
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is the momentum translation operator e™ \P) = 

ik)) 
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exp (— zKCOsX) = ^ (— z)"* J,„(K)e*' 



where we introduced the quantity 

K 



(4) 



Therefore, if one starts from an initial state of well- 
defined momentum |Po)- only states of the form 
\Po + mk) (m integer) will appear in the dynamics. We 
can write any arbitrary momentum Pq in the form 

Po = (m + I3)k 

with TO integer and (3 € [—1/2, 1/2) (i.e. /3 is in the "first 
Brillouin zone"). This introduces the quasimomentum kf3 
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which, as just demonstrated, is a constant of motion. The 
particle wavefmiction ^^{X) = {X\ ijj) can be written as 



1 



ij{k)dk, 



where k = P/k, and can thus be decomposed in: 

hl/2 



1 /• + 1/2 +°° 



,i{m+(3)X 



-1/2 



-1/2 



(5) 



where we introduced the quasimomentum component 

MX)-- 



i(m+/3)X 



(6) 



with 



The function 



(8) 



is of a 27r-periodic function. Eq. ([S]) is a direct manifes- 
tation of the Bloch theorem: the particle is described by 
a Bloch wave, il;p{X), which is the product of a plane- 
wave of well-defined quasimomentum /3 and a periodic 
function up(X). 

In most experimental realizations of the kicked rotor 
with laser-coled atoms the initial velocity distribution is 
larger than the Brillouin zone of the system and all quasi- 
momenta are present; it may thus be necessary to aver- 
age observable quantities over the quasimomentum. It is 
shown in the App. [A] that these averages values can be 
simply expressed in terms of the BW decomposition: 
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-1/2 
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(9a) 
(9b) 



(j^ where the subscript /3 indicates an average over "0^. 



In order to unveil the physical origin of quantum res- 
onances, let us take, for simplicity, an initial state of 
quasimomentum (3 = and apply to it the free evolu- 
tion exp(— i^P^/2): 
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ijQ{X,t^\-) = 



1 .r-^ ~ f ^'rn'^\ 

-== ipo ('71, t = 0) exp I —ik—^ j exp (imX) . 



(10) 



If we set in the above expression 




with £ arbitrary positive integer, we see that the argu- 
ment of the exponential phase factor in Eq. (fTU]) is al- 
ways an integer multiple of 27r. This means that the free 
evolution over a period leaves the wavepacket invariant. 
Fig. [1] shows an initial B W that is spatially localized [111 
in the first Brillouin zone and which evolves freely with 
time. The BW first becomes completely delocalized but, 
as seen above, it "focalizes" back to its initial shape just 
before the kick! Our using of an optics language is more 
than fortuitous: a monochromatic beam diffracted by a 
grating also reforms after a certain propagation length, 
an effect known in optics as the "Talbot effect". The 
quantum resonance is the atom-optics analog of the op- 
tical Talbot effect 

The only possible evolution comes from the kicks, 
whose effect is to add new momentum components ("side- 
bands") separated by a multiple of ^. As the effect of 



the kick is diagonal in position representation; it adds a 
position-dependent phase to the wavefunction: 

'ipQ{X,t = 1) = exp(-iKCOsX)Vo(^,^ = 0). 

As the wavepacket has always the same shape when the 
kick is applied, it acquires the same phase from each 
kick, and the kick effects add coherently. The quantum 
resonance condition is thus a constructive interference ef- 
fect (in contrast with dynamical localization for instance, 
which results from destructive interferences that "freezes" 
the wavefunction evolution). 

The above discussion puts into evidence an analogy 
(despite its very different nature) between quantum and 
classical resonances: in both cases, for particular values 
of the parameters, the dynamics is such that the particle 
(the wavepacket in the quantum case) is always kicked at 
the same position and the effect of the kicks adds (con- 
structively interferes) to produce a linear increase of the 
momentum. An important difference, however, is that in 
the classical case the effect is related to the intensity of 
the kicks - the resonance condition thus depends on the 
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(a) 



(b) 



(c) 



(d) 
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Figure 1: Free evolution of a localized wavepacket (/3 = 
0)between two successive kicks for k = 47r. (a) Initial 
wavepacket centered at X = n/2 shown in the interval 
X G [— tt.tt], (b) wavepacket at t = 0.115, showing complete 
delocalization, (c) wavepacket ai t = 0.25, showing two repli- 
cas of the initial shape, (d) "reconstructed" wavepacket at 
t = 1~ , identical to the initial one. 



pie" resonances" corresponding to ^ = 2nl {£ positive in- 
teger) are analyzed, in which the initial wavepacket "refo- 
cuses" in a single replica of the initial wavepacket. In sec- 
tion llVl we focus on the so-called "high-order resonances" 
(k = A'Kr/s, r, s integers), where the free evolution of an 
initial wavepacket generates various spatially-separated 
replicas of the initial wavepacket. 



III. ANALYSIS OF "SIMPLE" QUANTUM 
RESONANCES 

Quantum resonances obeying the condition k ~ 2t:£ 
(with £ a positive integer) arc named "simple" quantum 
resonances (SQR). With this form for k, it is shown that 
the shape of BW remains invariant in the free propa- 
gation step between two kicks and that the effects of 
the kicks may lead to a ballistic (linear) growth of the 
momentum; the kicking period is a multiple of the half- 
Talbot time defined in classical optics, which is the con- 
dition for the optical (integer) Talbot effect [s^, H^l • 



parameter K - whereas in the quantum case it is related 
to the constructive accumulation of quantum phases and 
the resonance condition thus depends on the value of k. 

In the following sections we shall develop an approach 
to calculate and interpret the behaviors of the system for 
different types of QRs. In section Hill the so-called "sim- 



A. Time-evolution 

Let us consider, at some (integer) time {t — 1) a, state 
^pj3{X,t—l) corresponding to the general form of Eq. 
and apply to it the free-evolution operator with k = 2Tr£ 



expf-i^^M^^t'^) = -^^#(m,i-l)exp(-z^£(m + /3)2)e^("+'5)^ 

= exp (— i7ri!/3^) ipp{m, t — 1) exp (— zTr^m^) exp {—i2n£m(3) 



gi(m+/3)X 



As Im^ and £m have the same parity, the first exponential factor under the sum is equal to exp {—iiTim) ; it can thus 
be combined with the second term, yielding exp {—imk {[3 + 1/2)). We obtain 



2 - 

i(m+/3)X 



1 / kp^\ ~ 

ipl3{X,t~) = exp ( ) V'/3("^7 ^ - 1) exp {-imk/S' 



with f3' = f3+l/2. This last expression can be rewritten, 
using Eq. ©, as 



V'^(X,i-) =cxp(^*M^±i^)v^0 iX-kf3\t-l). 

Applying now the kick operator (which is diagonal in the 
X representation), we obtain a recurrence relation linking 



the wavepackets at times t and (t — 1) 



Mx,t) 



(12) 

The above result shows that in the conditions of a 
simple resonance \il}{X,t)\^ — |?/'(X — w, < — 1)|^ = 
\'4'{X — ut, 0)1^: the square-modulus of the wavcfunction 
remains invariant immediately before the kick, except for 
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(a) 



(b) 




(c) 



(d) 
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and 



il)p{X-n,t= 1) = e*' 



which means that the phase added by the kicks simply 
cancels after two kicks. An analogous cancellation hap- 
pens for any other value of the velocity v obeying Eq. (jl6p 
with q ^ 1: the force acting on the particle averages to 
zero after tr kicks, and the motion of the wavepacket is a 
simple oscillation of period tj.. This behavior, that can- 
not lead to a ballistic behavior, is called "anti-resonance". 
If 1) = 2-Kp [or g = 1 in Eq. corresponding to 



Figure 2: Wavepacket plotted at different integer times, for 
k = An and an irrational quasimomentum (/3 = -^^ 0.063). 
(a) Initial wavepacket (IV'/jI^) centered at X = 0, (b) t = 1, 
(c) t = 2, (d) t = 5. The shape of the potential is represented 
in dotted lines. The wavepacket interacts with the potential 
at different positions, and the average effect tends to zero. 



the evolution of the BW is given by 



a drift with a "velocity" [4l| 



(13) 



Moreover, the wavefunction acquires a position- 
dependent phase due to the kick. 

By iterating Eq. down to i = 0, we can express 
ipf3{x, t) in terms of the initial wavefunction 



^Pl3{X, t) = exp t)) (X - vt, 0) , 

where the accumulated phase is given by [i^ l 



E 

s=0 



COS {X — vs) . 



A remarkable property is obtained if 

P 



2Tr- 



q 



(14) 



(15) 



(16) 



with p,q integer, e.g. for rational values of the quasimo- 
mentum. Then, the phase defined in Eq. (fT5|) takes the 
values ^{X, t = q) = ii q ^ 1 and <i>(X, t ^ q) ^ cosX 
ii q ~ 1. Let us first consider the case q ^ 1. After a 
recurrence time tr = q the BW is given by: 

^I3{X, tr)=ijfi{X~ 2ttp, 0) = e'^^P^MX, 0), 

and the particle comes back to its initial state after tr 
kicks, leading to a periodic evolution. To show the phys- 
ical origin of this periodicity, let us take the simple case 
V = TT, corresponding to = 2. From Eqs. ((M]) and 
(fT5)) . the BW evolution over two successive kicks is 



MX,t^2) 



-IK COS X 



ijjp (X~7r,t= 1) 



^f,{X,U = 1) 



-iK. COS X 



ibp (X - 2np, 0) 



g-»2,rp/3g-«cosJf^^(^^ 0). 



The wavepacket exactly recovers its initial shape after 
each kick and the particle is subjected to an identical 
potential at each kick. In contrast with to the former 
case, the kicks will add in a coherent way, which, as we 
shall demonstrate in Sec. IIIIBI causes a linear increase 
of the average momentum, or ballistic behavior. This is 
the quantum resonance. 

If Eq. is not fulfilled, our picture still allows to 
guess the general shape of the asymptotic evolution for 
large numbers of kicks (i ^ 1). If w/27r is not a rational 
number, the phase ^{X,t) will tend to as t ^ +oo. 
Once again, the contribution of the successive kicks will 
average to zero and the wavepacket just "drifts" (for stro- 
boscopic times) while keeping its initial shape (see Fig. [2]). 



B. Averages 

Let us now focus on the time evolution of average val- 
ues. From Eq. p2p we can easily obtain a recurrence 
relation for the BW average position 



(A) {t)^{X) {t-l) + v. 



(17) 



A recurrence relation for the average momentum is ob- 
tained using Eq. p^ : 



(P)^(t) = {P)^it-1)+K / dXsmX\i,,3{X -v,t-l)\' 



{P)g {t-l)+ K / dX sin(A + vt) |V'/3(A, t 



(18) 

This expression can be iterated down to t = 0, leading 
to: 
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(P),3 (<) = {P)p (0) + K ( ^ / dX sm(X - 



nv)\MX,0)\' 



(19) 



The role of interferences is clearly seen in this last result: 
indeed, if w = k{P +5) = Stt [mod 2tt], the momentum 
increases linearly with the (stroboscopic) time i, i.e 

{P)p{t)^{P)^j{0) + Dt (20) 

where the slope 

D^K'f dxsina;|V'/3(x,0)|^ (21) 

clearly appears as a force averaged over the initial spatial 
distribution. The observed momentum transport shows 



that well-controlled diffusion can be obtained by a suit- 
able choice of initial conditions [i.e of ip^ix, 0)]. 

For V = k{(3 + i) = TT [mod 27r]. destructive interfer- 
ence occurs. This can be seen from Eq. ([T^ : (-P)^ (t) = 
(P)^ (0) [t even) or (P)^ (t) = (P)^ (0) - D {t odd). In 
the general case {v ^ 27r[mod 27r]), the momentum is 
frozen around its initial value. 

Quantum resonance effects can also be analyzed 
through the temporal evolutions of the average kinetic 
energy. Using Eq. p^ . one obtains: 



(E)p{t) = {E)p{t = 0) + ^ y" dX sin(X + nv)^ li^piX, t = 0)\^ + K J dX sin(X + nv)j J{X, t = 0) 

(22) 

where we introduced the current 

J(X, t) = {i,p{X, t)dxi'}AX, t) - c.c) . 

The constructive interference case f = 27r [mod 27r] then leads to an average kinetic energy increasing quadratically 
with time: 

(E)^{t) = (E)^{t = 0) + ^Kh^ f dX sin^ X \'iPp{X,t = 0)\^ +Kt f dX smXJ{Xx,t = 0) (23) 



which is the quantum-mechanical analog of the ballistic 
motion observed - in different conditions - for a classical 
resonance. The ballistic growth is seen to be propor- 
tional to the quantum average of the square of the force 
sin'^X. 



C. Map 

Inspection of Eqs. (fTT]) and Eq. p8|) suggests that 
the dynamics of position and momentum averages of a 
Bloch-wave can be described by a map. In the follow- 
ing, we assume that the initial BW is sharply localized 
around its mean position {X)^ {t = 0) = Xq. Noting 
Pt = {P)p {t) and Xt = {X)fj {t) we have from Eq. (HH), 
Pt = Pt-i + K sin(Xo + vt) (with the normalization con- 



dition /^^ dX \tl;p{X,t^ 0)f = 1). We then obtain 

Pt = Pt.i + KsinXt (24a) 
Xt = Xt-i+v, (24b) 

which evokes the classical map (see Eq. [2]), with the 
important difference that the position in Eq. (|24b[) does 
not depend on Pt but solely on the drift velocity v. This 
produces after t kicks 

Xt = Xo + vt 

Pt = Pt-i+ Ksm{Xo + vt). 

Iterating in turn this momentum equation produces 

Pt = Po + i^ELi sin(Xo -f nv). 
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Ballisticity is found if w = 27r (mod 27r): 
Pt = Po + tK sill Xo = Po + Dt 

with 

D^A'sinXo. (26) 

For a localized packet, ballisticity emerges, as in the clas- 
sical case, if the particle is always kicked at the same po- 
sition. In the classical case, this is possible \i K = 2ttp 
(p integer) and Xq = ±7r/2. In the quantum case, the 
resonance condition depends on the quantum parameters 
k and P (through v), but the kick intensity K and the 
initial position Xq determine only the growth rate of the 
average momentum. 

Recalling that t is an integer (stroboscopic) variable 
counting the kicks, one can find also a whole class of 
periodic behaviors. An example is ^ = 27r and /3 = 
{v = tt), with a recursion time tr = 2, which has al- 
ready been analyzed in sect. IIII Al More generally, for 
rational quasimomentum values f3 = p/q {p,q integer), 
after a number of kicks tr-, the average momentum and 
kinetic energy come back to their initial values. In such 
cases the dynamics is periodic, due to the effect of kick 
compensation discussed above. 



with the situation found in simple resonances, after a 
free propagation the initial packet docs not reconstruct 
in an identical packet, but forms two or more replicas 
of the original one. The action of the kick in these sub- 
packets generates different quantum phases and produces 
quantum interference effects during the subsequent free 
propagation. This makes high-order quantum resonances 
fundamentally different from, and more complex than, 
simple ones. 



High-order quantum resonances (HQRs) correspond to 
a dimensionless Planck's constant of the form k = Anr/s, 
with r and s > 2 integers. A calculation similar to that 
leading to Eq. (112^ shows that after the free propaga- 
tion the initial packet refocalizes into s uniformly spaced 
subpackcts if s is odd, and s/2 subpackets if s is even. 
We shall consider here only the the simplest case k = tt. 
The method presented below can in principle be general- 
ized to more complicated cases, but the algebra involved 
quickly gets very cumbersome. As in the preceding sec- 
tion, our approach gives a simple picture of the physical 
mechanism of HQRs. 



IV. THE k = TT HIGH-ORDER RESONANCE 

"High-order" QRs are the quantum-mechanical analogs 
of the fractional optical Talbot effect In contrast 



Using the general expression for BW at time {t — 1), 
Eq. ([6]), and applying the free-evolution operator pro- 
duces in the case k = tt: 



MX,n = ^^V/3(n,t-l)exp(-z|(n + /3)')e^("+«^ 



cxp — I 2^ t - 1) exp — I exp [-iTTuf]) exp {i{n + P)X) . (27) 



We show in the Appendix |B] that the above expression can be written as 

-iK cos X 



V2 



[e-*T V/3 {X -kf3,t~l) + e'Te'l^''ijfi (A - S/3 - tt, i - 1) 



(28) 



or 



ipp{X + wt, t) 



-iK<p{X,t) |- 



V2 



e {X + w{t-l),t-l) + e'''^^e"^''^pf3 [X + w{t - 1) - TT,t - 1) (29) 
I 



where w is the packet (stroboscopic) drift velocity defined and we introduced the "local" phase [43] 

by 

(f>{x, t) = K cos{X + wt) 



(31) 



w = kp 



Eq. ([^5]) shows that at any time t, ipi3{X + wt, t) is the 
(30) superposition of two subpackets having the same shape 
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as the initial BW, centered at X = and X = w. Each 
of these two subpackets is multiphed by a phase factor 
which is the sum of the accumulated phases, producing a 
complex interference pattern. The principle of the follow- 
ing calculation is to keep track of the coefficient of each 
subpacket, as we know that the shape of the subpackets 
is fixed. We then write: 



iX + wt,t) = Ci (X, t)^f3 (X, 0) + C2 (X, t)ij(3 {X - 



-^,0), 
(32) 

where ci (X, t) and C2{X,t) are 27r-periodic complex am- 
plitudes |44l |. The above expression corresponds to a 
coupled two- level model where the "particle", initially in 
"level 1" (i.e ci = 1 at i = 0) is progressively "trans- 
ferred" to level 2" (c2 7^ 0) and then back again to "level 
1", performing a kind of Rabi oscillation. 
Let us define the state vector 



ciiX,t) 

C2(X-7r,<) 



It is shown in App. IB 21 that these amplitudes obey a 
matrix recurrence relation : 

Ct = MfCt-i 

where Mt is a matrix depending on time and space having 
the form 



where 



~ 1 



V2 V *e'^e 



(33) 



with (l){X + wt) given by Eq. ([?!]) . The matrix Alt can 
be recast as 



M, 



1 



V2\ 



Q 



1 ie-'f^'' 



1 



The rightmost matrix in the above product stands for 
the free propagation that induces a coupling between the 
two subpackets; it is thus responsible of the interference 
effects. The leftmost one represents the effect of the kick 
and is obviously diagonal in x-representation. 

Analytical results can be obtained in the case /3 = 0. 
As in this case w = 0, the matrix Mt is time- independent. 
It is then easy to write Ct as a function of the initial 
condition 



Co 



Ct = [A/t]* Co 



The eigenvalues of Mt are: 



Mt 



Co 



A = exp(±i6) 
where the phase 9 depends on X and is given by 



cos Q 



cost 



n/2 

cos(kcosX) 

(note that 7r/4 < < 37r/4). If P is the diagonalizing 
matrix, one can write (see Sec. IB 31 in App. |B]) 



^0 g-.te ) P-'co. (34) 



After a straightforward calculation, the amplitudes are 
found to be 



, ^, i\/2 sin0 . , 

cos(te) - — ^ sin iO) 

2 sni B 



(35) 



and 



C2(X-^,i) = ^-^e-**7 sin(te). (36) 



From the above result we can calculate the average mo- 
mentum, but the algebra involved is cumbersome (see 
App. P: 



{P),,{t)^{P)^{t-l) + K / dXs\xiX\ijp(X,t)\\ 



In the general case, it is difficult to analyze the be- 
havior of the amplitudes ci {X, t) and C2 {X, t) since they 
depend implicitly on X in a quite complicated way. An- 
alytic results can however be obtained in the limit which 
the initial BW is formed of spatially narrow wavepack- 
cts.Lct us assume that ■0^(X, t = 0) is well localized at 
some position Xq in the interval [— tt, tt] (i.e its width 
AX ^ 27r) and express the average momentum in terms 
of the coefficients ci and C2 



{P)p{t) - (P)^(t-1)+X / dxs\ii{X + wt)(\ci{X,tf -\c2{X~^,tf)\^p{XM 



(37) 
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This expression shows that if the two subpackets have 
the same weight, the momentum shift per kick is zero: 
the subpackets are locahzed around positions Xq and 
Xq + TT, and subjected to opposite forces +K sm{Xo)/2 
and —Ksm{Xo)/2, an effect that is characteristic of the 
high-order resonances and obviously does not exist for 
simple resonances. This expression is valid if the two 



subpackets are well separated so that they do not inter- 
fere significantly. 

For a wavepacket that is strongly localized at Xq, the 
amplitudes in the decomposition Eq. (|37p can be evalu- 
ated at Xq and depend only on time, while the phases (fi 
and Q take constant values {(f> = kcosXq). Hence, the 
average momentum evolution is 



(P)^(t) = {P)^{t-l) + Ksm{Xo + wt) (l-2|c2(Xo-^,i)|') 
which, in our simple case /? = 0, takes the explicit form 



(38) 



(O,=o it) = iP)p=o it-l) + KsmXo[l- 

^ sm B / 



Note that the expression inside parenthesis in the above expression is characteristic of the diffraction on a grating, 
which puts into evidence the "wavelike" nature of the dynamics. Iterating down to t = produces 



{P)p=o it)^{P)p=o (t = 0) + XsinXo 



sm 



1 + sin-" 

r 



t + 



1 



1 + sin^ 



sin[(2t+l)e] 1 
2 sine 2 




-71 -7t/2 7l/2 71 



Figure 3: Momentum slope D [Eq. (|39|l ] as a function of Xq 
for 3 values of k: 0.5 (full line), 1.0 (dots) and 2.0 (dashed). 



This expression shows that average momentum evolution 
is, for arbitrary cp and &, a mix of two qualitatively differ- 
ent behaviors: ballisticity, corresponding to the first term 
in the brackets, and oscillation, described by the second 
term. This contrasts with SQRs, where the dynamics is 
either ballistic or oscillatory. 

For long times and 7^ 0, the ballistic term dominates 
and one gets 

{P)b=o (i) = (^>fco it^O) + Dt 



with 

D = KsmXo( "^""'t V (39) 

being the rate of change of the mean momentum. The 
behavior of D as a function of Xq is displayed in Fig. [3] 
for different values of k. It is interesting to compare 
Eq. ([55)) with its counterpart for SQRs, which is given 
by Eq. (P5|) : note in particular that in the present case 
D = if Xq = 7r/2, whereas for SQR, the maximum of 
D occurs when the force is maximum (i.e at Xq = 7r/2). 

The coefficient D depends periodically on the kick in- 
tensity K via (j) = HL cos Xq, as shown in Fig. 3] In con- 
trast to the simple resonance case, increasing the kick 
force does not necessarily increase the diffusion, an ef- 
fect that persists if the momentum is averaged over the 
quasimomentum distribution. 

Analogous results describing the ballistic behavior of 
an initial state which is an eigenvector of the momentum 
have been obtained via a quite different approach [i^ . 
Our localized-packet approach provides a clearer picture 
of the underlying physics. 

For /3 7^ 0, the occurrence of ballisticity depends on the 
periodicity of Mt and on the initial conditions. More pre- 
cisely, ballisticity will emerge if the relation Mt+t^ = Mt 
is fulfilled for some (integer) recurrence time tr- This 
happens for any rational value of quasimomentum. For 
example, if /3 = 1/2, Mt [see Eq. ([55)1 ] has a period of 4 
kicks: Mt+4 = Mj. One can apply the above reasoning to 
M = M4M3M2M1. To illustrate the resulting behavior, 
Fig.Oshows the time-evolution of average momentum ob- 
tained by direct integration of the Schrodinger equation. 
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One observes essentially the same kind of behavior. 



For a general rational quasimomentum f3 = p/q, the 
ballistic diffusion rate is roughly proportional to 1/ q. Ir- 
rational quasimomenta, that may be consider as the limit 
q +00, do not produce ballistic behavior. 



1 2 3 4 5 



Figure 4: Momentum slope D [Eq. (|39[l] as a function of k for 
two initial positions Xo=7r/4 (dots) and Xq = tt/S (full line). 
One sees that increasing k (that is, the kick intensity) does 
not necessarily increase the momentum slope. 




-1 I • • • • ' • ' 1 

5 10 15 20 25 30 35 40 

t (kicks) 



Figure 5: Evolution of the average momentum for /3 = 1/2 
obtained by numerical integration of the Schrodinger equa- 
tion. The initial initial wavefunction has the same Gaussian 
shape as in Fig. [T] and is centered at the positions Xq = 7r/4 
(dark line), which displays a dominant oscillatory behavior, 
or Xq = 7r/8 (light line), displaying dominant ballisticity. 



V. AVERAGING OVER QUASIMOMENTUM 



In experiments performed with laser-cooled atoms (not 
with ultracoldsdoms - BECs) the initial momentum dis- 
tribution is larger than the Brillouin zone unless veloc- 
ity selection is performed. In position space, this cor- 
responds to an initial wavepacket which is localized on 
a single potential well, around a position Xq. In this 
simple case, and for k = we are able to give an an- 
alytical expression for any value of of the momentum 
and kinetic energy averaged on quasimomentum /3. 



Starting from Eq. ([5a|) : 



{P)(t) ^ {P){t^Q) + Ky^ dp dXsin{X + nv)\i:p{X,Q)\^ . 

Assuming that V'/3(-'^7 0) is /3— independent (V'/3(^, 0) = ip{X)ai\d ^'^^ dx \ip[X)\^ = 1 ), and recalling that v — 
2tt£{P + 1/2), we easily see that averaging on (3 lead to: 

(P) it) = {P) it = 0). 

A similar reasoning can be applied for kinetic energy. From Eqs. (|9bp and (|22p. we have: 

{E){t) = {E){t = Q) + — dp dX sm{X + nv)] \(p{X)f + kJ^ dP dx sm{X + nv) J {X , t 

2 J -1/2 Vn=l / „=l^-l/2 
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When wc integrate over /3, the term of the last term cancels out, and the only contributions come from 

J ' ^X] ^"^(^ + " / ^ sin2(X + -Kin + 27r^??,/3)^ = ^ 

I 



We finally get the diffusive behavior 

{Em = {E)[t = ^) + ^. (40) 

This result evokes the classical kicked rotor in chaotic 
regime, whose kinetic energy grows linearly with time 
with the same rate. This shows that the ballistic behav- 
ior, which corresponds to a "null- measure ensemble" of 
rational quasimomenta, is very hard to detect by mea- 
suring quantities averaged over quasimomentum. 

Experimentally, the optimum situation for observing 
QRs is to perform a quasimomentum sclcctioix, either by 
using stimulated Raman transitions ,29. , 46 . .471 or by us- 
ing a Bose- Einstein condensate [13, lUHlj. However, it 
is possible to detect QRs with atoms issued of a magneto- 
optical trap if one can measure the full momentum dis- 
tribution with enough precision to see the ballistic parts 
of the wavefunction separating out of the diffusive part 
for long enough times, as experimentally evidenced by 
d'Arcy et al. [2l|. 

VI. CONCLUSION 

In the present work, we presented a description of 
quantum resonances of the kicked rotor in position space, 
both for simple and for high-order quantum resonances. 
We have shown that the dynamics can be understood by 
considering that the spatial wavepacket comes back to its 



in ital form after a finite number of kicks, according to the 
(rational) value of the quasimomentum. For a localized 
wavepacket (or an finite ensemble of localized wavepack- 
ets, in the case of HQRs), one can interpret the dynam- 
ical behavior in terms of the action of a finite number 
of successive kicks. This picture, inspired of the atom- 
optics analog of the Talbot effect, proves very useful both 
as it providing a intuitive understanding of the underly- 
ing physics ans as it leads to analytical developments for 
experimentally-relevant quantities. 



Appendix A: AVERAGE MOMENTUM 

The average momentum reads: 

(P> = ^V/ d/3(n + /3) V/3 (n) 

n -^-1/2 

and can be related to {P) ^ in the following way. One has 
from the definition of section [III 

V27r J_7r 

therefore: 



(n+/3)^/3 (n) = 



OA v27r 



# (X) , 



The first term on the RHS vanishes {up {X) = ipp (X) e '■^^ is 27r— periodic). One then has: 

ik 



ik 
'2^ 



d/3 



d 



— TT J — TT 



dXdX'e-^^^^'-'^'Hp [X') {X) ^ e 



Finally, using Y.n e-^<^-^"^ = Y.k '^(^ - X' - S/ctt) 



(p) = ~ik I dp / dxrp {X) —ijp (x) 



1/2 



1/2 



d/3 (P) 



where 



r dxrp {X) [Pi^p {X)] 

J — TT 
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For the kinetic energy, the same developments lead to: 
where 

{P^)p = r dX [P^p (X)]* [P4'0 {X)] 
(note that the same rules can be obtained for (P^) = / d(i{P^)^). 

Appendix B: BLOCK WAVE EVOLUTION FOR THE % = -k RESONANCE 

1. Bloch wave evolution 

The free-propagation factor exp(— i|-n^) in Eq. ([?71) has two different values according to the parity of n: 



Replacing 



exp (~*2''^^) ^ ^ in even) 
= —i (n odd). 



i;p{n,t)^^f dXe-"'''up{X,t)^ [ dXe-*("+«^#(X,t) 
v27r j_7r v27r J_7r 



in Eq. (gT]), one finds 

exp (-i7r/3V2 



2tt ^ 



one can separate even and odd terms ; this leads to: 



CXp_^^^V2) J\^,^^^^,^^_^^^-^fJiX-X') Y^^-^2.pP^-^2p(X-X')(^^ 



Using the relation e *27rp/3g i2p{x x ) _ (1/2) 6{X — X' — 7r/3 + n7r) and integrating with respect to X' , gives 
Eq. ([291) (note that only n = 0, 1 contributes for /3 > and n = 0, -1 for /3<0). 

2. "Two-level" system 

The coupled equation for the amplitudes ci 2{X,t) are obtained in the following way. Insertion of Eq. (|32p at time 
(t~l) 

ijp{X + w{t -l),t-l) = ci{X,t- l)^p{X, 0) + C2(X, t - l)i^p{X - ^, 0), (Bl) 

in Eq. ^ gives 

#(X + wi,t) = ^^^^ (e-'"/4[ci(X,t-l)V'/3(X,0) + C2(X,t-l)iA^(X-^,0)] + 

g,7r/4g»/3,r (-^ _ 7^, i - f3{X - TT, 0) + C2 (X - TT, t - l)V'/3(^ " StT, 0)]) 
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and can be put in a simpler form [using ^Jf}{X — 27r, 0) = e ^'^^^il)p{X^ 0)] 



V2 



( [e-'"/4ci(X, t - 1) + e}^/^e-'^''c2{X - tt, t - 1)] i^p{X, 0) + 



e*'^/4e*^'^ci(X-7r,t- l) + e-*'^/^C2(X,t- 1)] V/3(^-7r,0)) 



Comparing to Eq. ([32]) one obtains: 

C2(X,i) 



e~"/-*ci(X,t - 1) + e"/^e~'^''c2{X - 7r,t - 1) 
e^^/'^e'^'^ci (X - TT, t - 1) + e''''I^C2{X, t - 1) 



This last expression can be put into the form 



C2{X~^,t) 



V2 



gW4g»/3^cj(j^ - 27r, i - 1) + e-*'^'"'c2(X - tt, i - 1) 



i7r/4 



One can show that the amplitudes are periodic functions, i.e ci 2{X, t)~ ci^2{X+2tt, t) by combining Eq. (j32p combined 
with the equality tppiX — 27r,t)= e~'^'^^V'/3(-''^: This property leads to the matrix expression 



ci{X,t) __ 
C2{X + Tr,t)J ^ 



ci{X,t-l) 

C2{X + ^,t-\) 



3. The case /3 = 



The explicit expression for the amplitudes ci^2{X,t) is obtained for /3 = as follows. The diagonalization matrix 

P is formed with the eigenvectors of Mj. For A = e^^'^ they are given by (— ie"*"^, e"*"^ — \/2e^*®)-'" . P is then 
obtained as 



V2e-'Q 



In order to obtain the amplitudes ci and C2 at time one performs explicitly the development corresponding to 
Eq. fM]) . The algebra is simple, although rather long, and the final result is Eq. 



Appendix C: CALCULUS OF AVERAGE VALUES FOR THE £ = tt RESONANCE 

Starting from Eq. (j28p , recursion relation for ?/)3 {X, t) , we easily obtain a recursion relation for its derivative 

^V/3(^,0 = ^KsinXe-'''™^^V/3(^,^) 

+ ^ (^e-^— V/3(X-^,t-l) + e'^e'^-— V;3(X-^;-7r,t-l)J . (CI) 

Using these expressions into for calculating the average momentum produces two terms; let us call them pi and p2- 
The first one is given by: 



pi ^ K I dXsY[YX\il),3{X,t)f 

J — TT 
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For the second term p2, on obtains 

P2 ^ dX {e'^il;j;i{X -w,t^l) + e~'^e-"'^''i;*p{X -x-Tr,t-l)) 

X (^e-''^ ^i^p (X-w,t^\) + e'T e'/^. - u; - ^, f - 1)^ 

+ \ j^^ dX [e^f^-r^ {X, t 1) (X - ^, i - 1) - e-'^V; (X - TT, t - 1) (X, t - 1)^ 

[in the last equahty we replaced (X — 'w)hy X, as the integration is over one period, the integration limits can be kept 
the same]. One recognizes (P) ^ — 1) in the first term while the second term cancels out. The recursion relation for 
the momentum thus reads 

{P)p{t)^{P)p{t^l) + K [ dXsinX\^l^p{X,tf . (C2) 
We now use the decomposition of ■0/3 proposed in Eq. (I32p . Eq. (jC2[) then transforms into: 



(P)^(i) = (P)^(t- l) + iA'cos(/37r) / dXsin{X + wt){c*{X,t)il:}{X,Q) + cl{X,t)TP}{X -TT,Q)) 

X {ci{X, t)i:p{X, 0) + C2{X, t)MX - TT, 0)) . 

If the initial wavefunction has a narrow distribution centered around position Xq, all terms involving overlaps of 
functions V/3 i^, 0) and {X + n) tend to zero [11] . This expression therefore simplifies to 



(P)^(i) = {P)p{t-1)+Kj dX siniX + wt) [\ciiX,t)f -\c2iX-TT,t)\'')\ijpiX,0)\'' (C3) 
Finally use the fact that ijjp{X, 0) is much narrower then any other factor, and assuming J^^ dX \^f-j{X, 0)|^ = 1 [ioj 



=cives 



(P)^j(<) = (P)^,(t-l) + Xsin(Xo + ^«t) (l-2|c2(Xo-^,t)|') 



where we used the normalization |ci(X, t)| + |c2(X — "'j'^)! = 1- 
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